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Abstract 

We introduce a new matrix model that describes Causal Dynamical Triangulations (CDT) in two dimensions. 
In order to do so, we introduce a new, simpler definition of 2D CDT and show it to be equivalent to the 
old one. The model makes use of ideas from dually weighted matrix models, combined with multi-matrix 
models, and can be studied by the method of character expansion. 



1 Introduction 

It has long been hoped that a theory of quantum gravity could be defined as the continuum limit of a 
discretised gravitational path integral. Using the Dynamical Triangulations (DTs) discretisation, the histo- 
ries to be summed over are "triangulations" built by gluing together equilateral simplices in every possible 
configuration, weighted appropriately for gravity [1]. With this simple definition, many results of physical 
interest can be obtained by analytical means. In 2D, matrix model techniques |2j have been particularly 
useful in this regard, allowing the solution of many models of 2D quantum geometry and matter. 

In 4D the original DT model does not seem to give physically acceptable results. However, adding a 
physically motivated restriction on the space of configurations leads to the better-behaved Causal Dynamical 
Triangulation (CDT) models [3]. Here, all the observables that have been calculated are consistent with the 
emergence of an extended 4D geometry at large scales [H [5]. The ability to calculate genuine generally 
covariant observables from a theory of quantum gravity, with good results, is a unique and important 
achievement for this approach. 

However, not all of the analytical tools that were so useful in the case of DTs have been successfully 
carried over to the case of CDTs. In the new model, the space of configurations is limited by demanding that: 
(a) there exists a global time foliation; (6) that spatial topology change is disallowed. At first glance, these 
"causal" conditions seem to be of an irreducibly global nature. Many of the analytical techniques used for 
DTs are not easily compatible with the imposition of such a global condition. Other techniques have found 
a role in 2D CDTs, which indeed has been solved, in the pure gravity case, by several different means [Si [7]. 
Unfortunately none of these methods has proved useful in the presence of matter (where the only analytical 
results come from high-temperature expansion [8]), or in higher dimensions, with the only exception of the 
heap of pieces representation [9] of a reduced model of 3D CDT [10] . 

In the original 2D DT model, the partition function can be re-expressed as a matrix integral, which in 
many cases can be explicitly computed. The techniques have been generalised in many ways {e.g. adding 
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matter, different face shapes, and non-standard weights). However, it has until now been thought impossible 
to impose the CDT conditions directly on a matrix model, and applications of matrix models to CDTs have 
followed other ideas [TTJ [T^ [T3j . Restrictions on the triangulations generated by a matrix model typically 
depend only on the properties of local structures such as vertices and faces. It is not obvious how restrictions 
of this type could be used to introduce global foliations and ban topology change. 

In this paper, we exhibit a novel matrix model that imposes the causal conditions, and generates CDT 
configurations. As in the original DT case, the correspondence between the matrix model and CDTs applies 
directly at the discrete level. The fact that this is possible is interesting in itself, opening up new possibilities 
for the study of matrix models and CDTs in general. Using this matrix model, we can address problems 
in the CDT model, such as finding the critical point. The results of these calculations will be presented 
elsewhere. 

The observation that makes a matrix model formulation of CDTs possible is that global restrictions, 
such as topological restrictions, can sometimes be brought about as consequences of local restrictions, by 
what might be called "rigidity" . For instance, if a finite tessellation is restricted to have only four-sided 
faces, and only 4-valent vertices, it is easy to see that the only allowed topology is the torus, since the 
tesselation must be a regular lattice. This rigidity has been noted in the context of matrix models before, 
where "dually weighted models" have been developed that allow non-trivial weights that depend on both 
the valency of vertices and the number of sides of faces [Mj [151 [K] ■ The CDT restriction can be expressed 
in a similar way, as a "partial rigidity". In a 2D CDT, each vertex in the triangulation is incident on exactly 
two spacelike edges, that connect it to its neighbours in its ID spatial slice. This property, along with the 
spacelike/timelike colouring of the edges bounding the triangles, completely characterises 2D CDTs, as we 
prove below in section [2l 

Expressed in this new way, the CDT condition can be imposed on a matrix model (see section [S]) . The 
model generalises the techniques of dually weighted matrix models to impose the first restriction, and it is 
a 2-matrix model so that it incorporates the spacelike/timelike labeling. We envisage a new application of 
the character expansion methods mentioned above as a possible way to solve the matrix model; some steps 
towards the solution are presented in section [3. II 

2 A new way of defining 2D CDTs 

The following definition of 2D CDTs follows from the standard one [6] by replacing the terminology of 
triangulations with that of the dual fatgraph. Vertices in the triangulation are dual to faces bounded by 
loops (or "boundary-components") in the fatgraph as usual; now, spacelike edges in the triangulation are 
dual to timelike edges in the dual and vice-versa. Below we will use the colour labels A and B for spacelike 
and timelike (the same names that we later give the corresponding matrices in the matrix model). In a 
CDT, triangles, which had two timelike edges and one spacelike edge, become "AAB vertices" incident on 
two spacelike and one timelike edge in the dual. 

A CDT fatgraph is a finite connected fatgraph obeying the following conditions. (1) All vertices are 
AAB. (2) There are disjoint sets of faces, which we call "strips", each of which is connected, and such that 
its members are glued on their timelike edges only to other members of the set (these are dual to cycles of 
spacelike edges in the triangulation). Each face in the strips in the dual fatgraph has exactly two timelike 
edges glued to the adjacent faces in the same strip. (3) It is possible to label these strips sq, si, st-i, 
such that the following holds: each face li (with labels natural to the cyclic order of the strip) on strip st is 
glued by spacelike edges to faces on strip st+i (this addition is mod T for toroidial topology, while the rule 
is ignored for strip st~i for a cylinder). These gluings exhaust all of the spacelike edges (apart from the 
initial and final strips in the case of a spherical CDT, which are then glued to an initial or final face bounded 
entirely by spacelike edges The original description of 2D CDTs in mentions other properties, but they 

^In this definition, there is nothing to stop self-gluings and double-gluings between vertices, so a spatial slice may have length 
2 or 1. Also we have included an initial and final face, dual to an initial or final vertex rather than a boundary, so that these 
CDTs are spherical in topology rather than cylindrical. 
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(b) ' — : 

Figure 1: (a): A triangulation with time slices, and with triangles with 2 timelike edges and one spacelike edge, but 
with spatial topology change. Timelike edges are shown in red, and the slashes indicate identified edges. The same 
triangulation is shown twice to illustrate both the time slices and the fact that, in time slice t, there is a branching 
vertex incident on more than 2 spacelike edges, (b): The dual fatgraph of the same triangulation. 

follow easily from those given here. 

This formulation is not suited to the application of matrix model techniques. For this reason, we introduce 
an alternative definition: a CDT fatgraph is a finite connected fatgraph with only AAB vertices, and such 
that every face has either 2 or timelike edges. Below, we will show that the new definition is equivalent to 
the original CDT definition. 

One of the crucial features of the definition of 2D CDTs is the lack of spatial topology change: each 
time slice consists of one connected strip. The new definition reveals that a more local restriction, that is, 
a restriction on the colourings of vertices and loops in the fatgraph, can be used to enforce this topological 
property. In order for there to be topology change, one of two things would have to be allowed. If all 
vertices are to remain AAB, the only way is to allow "discrete Morse points" (or "branching points") in the 
triangulation [6] . As shown in figure [U such a branching point has more than two spacelike edges incident 
on it, leading to a face in the dual with more than two timelike edges. Another way of introducing topology 
change would be to add branching points in the dual, so that one strip would be glued to more than one 
strip in the same time direction. This cannot happen without there being a vertex incident on more than 
2 spacelike edges; the fact that all vertices are AAB is therefore also crucial. The comments below indicate 
how this CDT gluing between strips follows from the new definition. 

It can immediately be seen that a fatgraph satisfying the original CDT definition also satisfies the new 
CDT definition given above. The converse is not so obvious; to show this, we must take a fatgraph F 
conforming to the new definition and show that it obeys each part of the (lengthier) original definition in 
turn. Part (1) oi the definition is also a part of the new definition. As for part (2), in such a fatgraph F 
the set of faces containing two timelike edges can be partitioned based on which faces glue on timelike edges 
(this is dual to partitioning vertices by connectivity on spacelike edges), which is a partitioning into strips 
as described above. This leaves part (3) of the original definition, which is less trivial. 

Let us consider the topology of a strip s in F, considered as a set of faces joined only by the gluings of 
their timelike edges. The boundary is then the set of all spacelike edges in the strip. Because every fatgraph 
corresponds to an orientable surface (see e.g. [TTl [2]). and each face is connected to only two others in the 
strip, the strip has the topology of a standard strip (as opposed to, say, a Mobius strip). The boundary of 
such a strip is composed of two cycles. This observation allows us to show that the fatgraph does obey the 
third part of the original CDT definition. 

First let us treat the case in which F contains no completely spacelike faces (i.e. faces containing no 
timelike edges). In this case, all spacelike edges are in the boundaries of strips. Each strip is glued to other 
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strips (or to itself) by its spacelike edges, which, as previously shown, are partitioned into two spacelike 
cycles per strip. If one such cycle was glued to spacelike edges in more than one other such cycle, three 
spacelike edges would somewhere meet at a vertex, contradicting our assumptions. Thus, each strip is glued 
to only two others. The collection of strips themselves, with the relation between them of gluing, must form 
a cycle (they cannot form multiple cycles or otherwise the graph would be disconnected, which is disallowed 
by assumption). We can now choose a labelling of strips sq, si, st-i consistent with this cyclic order on 
the strips. These gluings exhaust all of the spacelike edges. This completes the proof that these fatgraphs 
obey the original CDT definition - more specifically, the definition of toroidial CDTs. 

There remains the case in which T contains some completely spacelike faces. The boundary loop of a 
completely spacelike face is a spacelike cycle. Again, this must be glued to exactly one other spacelike cycle. 
It is easy to see that only 2 such faces with no timelike edges can be included in the fatgraph if it is to be 
connected. Hence we can repeat most of the argument from the former case. In this case, it is possible to 
label the strips Sq, Si, st-i, such that the following holds: faces in strip St are glued by spacelike edges 
to faces on strip St+i, apart from when t = T — 1. These gluings exhaust all of the spacelike edges, except 
those from strips sq and St_i, which also glue to two completely spacelike faces. Thus, these graphs fulfil 
the definition of spherical CDTs given above. 



3 The matrix model 

We are now going to define a two-matrix model that generates, via its perturbative expansion, the fatgraphs 
discussed in the previous section, and whose partition function Z will be related to that of the CDT model 
in the usual way, i.e. 

ZcDT — fim lii^ . (1) 
The two-matrix model is defined by its partition function 

Z= [ dAdi?e-^Tr[iA^+i(c-iB)=-<,A=B] ^ (2) 



with C an external matrix satisfying (in the large N limit) the condition 

Tr (C") = N52jn for m > 1. (3) 

It is not hard to see that such model generates graphs for which in each face there are always two (or zero) 
_B-lines, and as explained in the previous section this is enough to impose the "CDT condition" on the dual 
triangulation. From standard analysis we can read off the free {g — 0) propagators of the matrix model: 

{A,jAm)o = ^SuSkj , (4) 

{BijBki)o = -^CiiCkj , (5) 

{AjBki)o - , (6) 

from which we find that in the expansion in Feynmann graphs of ([2]) a face contributes a factor of Tr (C™) 
with m being the number of internal i3-lines, and hence we can control such number by imposing conditions 
on Tr (C™), as in ([3]). Note that this is a generalization of the standard case where, being C = 1, a face 
contributes just a factor N. On the other hand and just as usual each edge brings a factor and each 
vertex a factor gN. In both cases the factors N gather to give for each graph F a global A^x(r) -^^ere 
x(F) = F — E + V = 2 — 2h is its Euler characteristic {F, E and V are respectively the number of faces, 
edges, vertices of the graph and h the number of handles of the surface on which it can be drawn), so that 
the expansion in is effectively a topological expansion, as first shown by 't Hooft [TS] . 
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The model can be simplified by noting that we can perform the integration over B, because it is a simple 
Gaussian, and obtain (up to a multiplicative constant which we discard from now on in the expression of 
the partition function) 

^ = |dAe-^Tr[iA^-4(A^cn ^ (7) 

which looks similar to the standard one-matrix model which describes two-dimensional Euclidean quantum 
gravity and which was originally solved in [19j (note that this kind of reduction from two matrices to one 
matrix was considered in [20^ for the case C = 1 showing that a model with Lorentzian triangles without 
further restrictions would be in the same universality class as the standard Euclidean model) but with the 
important modification in the vertex which turns out to be crucial. One can understand the integration over 
the matrix B as the gluing of triangles along their spacelike edge, which can be done only in one way, and 
which gives rise (if there are no boundaries) to a model of squares with only timelike edges. The presence 
of the external matrix C introduces an anisotropy for the squares, while the condition ([3]) acts to rigidly 
transfer such anisotropy from local to global level, i.e. to the whole quadrangulation. 



3.1 Character expansion 

The standard reduction to eigenvalues is not directly available in our case because, as [A, C] ^ 0, the two 
matrices are not simultaneously diagonalizable, so we need some fancier trick to solve the model, like the 
character expansion |14[ I15[ I16j . Indeed in terms of reduction to eigenvalues our interaction term poses 
exactly the same problem as the ABAB interaction of [16pl . i.e. we have to perform the integral 



^^giVVTr(OA^ntc)^ _ (8) 



where A is the diagonal matrix of eigenvalues of A. Unlike the case of the Itzykson-Zuber-Harish-Chandra 
integral (which differs from our problem for the absence of the overall square in the exponent) there is no 
exact formula for this kind of integral. Following [16] we can instead perform the character expansion 



{h} ^ [UW-l^V- 

where the sum is over all integers hi such that hj^ > /iAr_i > ... > hi > and which have to be in equal 
number even (h^) and odd (/i°), X{h}{x) is the character of the group element x (the group is GL(N) in our 
case) in the representation labelled by the set of integers {ft-}, = '}2,hi — ^N{N — 1) is the total number 
of boxes in the corresponding Young tableau, and A(x) is the Vandcrmonde determinant. 

Now we can perform the integration over the unitary group using the character orthogonality relation 
J (ifi X|/i}(r2A^ri'l'A(7) — X{h}{^'^)X{h}{C) /d^i^y (where is the dimension of the representation given by 

d{h} = A(/i)/n^7^«0> to obtain 



{h} 



where c^hy is the following coefficient ([a;] denotes the integer part). 



^Just replace their A with our and their B with our C. Of course our is a completely different model, because of 
the interaction containing four powers of A rather than two, because C is an external matrix, and because we have no other 
self-interaction term; but what we are interested in here is the similar xyxy structure which we want to disentangle. 
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and I^f^y is the matrix integral 



^JdA X{,}(A^)e-^^*^ ^ )>o , (12) 

which we have rewritten in terms of the fuU matrix for esthetic reasons (the angular integration is now just 
an overall constant, and hence we have the freedom to move it in and out of the partition function). Using 
the formula 

originally given in ^15] , for the character of a matrix satisfying the condition ([3]), we can rewrite Z as 

{h} 

We have now succeeded in reducing the original integral in 0{N'^) variables to a problem in 0{N) variables, 
and this is a promising step for its solution. 

Unfortunately the integral is not of the Di Francesco-Itzykson type [13] , and an exact formula is not 
known. Another possible way to tackle the evaluation of the partition function l|10p could be a double saddle 
point method along the lines of [16] . but again the quadratic dependence on A in the character makes the 
existing results difficult to use, and we could not obtain sensible saddle point equations. 

We could of course simply rely on the equivalence of the model to usual CDT, as proven in Sec. [21 and 
hence conclude that the model has already been solved. Although that would be a completely legitimate 
attitude, we find it instructive to try and take some further steps within the character expansion method, 
especially in view of the fact that, if successful, this method might be extended to solve other CDT models 
that have not been solved by other means. 

In order to proceed let us then look for a moment at the partition function Zi for the usual matrix model 
with quartic vertex, which is obtained by substituting C = 1 in Q, and which was first solved in J19!. If in 
([TUj) we substitute C = 1, wc find 

2 

fl5l 

- l^'^^T' rr r^.li (^{'^}(^ ' 

which of course in this expansion presents the same difficulty as our new model. In order to find an expression 
for (X{h}(^^))o we compare ([TS]) to the expansion that was derived for the same model in [15], 

E (^)*ri^™^n('.f'-'.f)(M"-'.f). (-) 

{/l(0),?l(l),?l(2),?i(3)} ^ ' £=0 Hi (-"^ )' »J 

where the integers h have factored into 4 groups of integers /i^'-' with e = 0, 1, 2, 3 denoting the congruence 
modulo 4. Matching powers of g among (jisp and (|16p . we can conjecture 

(X{m(^'))o -^Vi A'{2h^^^) l[{2h['^)U , (17) 

£=0 i 

as the extension of the Di Francesco-Itzykson formula we were looking for. We stress that this is by no 
means a proof, but just a conjecture that we would like to test against the known solution of CDT. 
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Substituting (|17p into our partition function (|14p we reduce it to a pure sum over representations, 



^ V2J 11 7 — Tz, — ^ n (z^?-/^?) ■ ^ ^ 

{h(0),h(i),h(2),h(3)} e^O ( n ( 4 " )! ) ^^'J^ 

We can now proceed to try and solve this model via a saddle point method for the shifted highest weights. 

Similarly to what is usually done in the case of the eigenvalues, in the large N limit it is useful to rescale 
hi — > hi/N, introduce the density distribution p{h) = for which we assume (following [HI [16]) that 

pih) = 1 , for < ft- < /ii , 

(19) 

< p{h) < 1 , for fti < ft < ^2 , 

for some hi and ft2 to be determined (note that the indices 1 and 2 have nothing to do with the index i 
which is omitted in the continuous notation) . This condition is simply the statement that in a typical Young 
tableau of N rows there will be a lower fraction of rows which are empty. Finally we define the resolvent 
function 

J Pih') h , f'- p{h') 



which as usual satisfies 

${h,) = lAiogA(ft) , for ft, e [0,ft2] , (21) 

where we used the notation -^(ft) = ^{H{h + ie) + H{h — ie). Assuming that, because of symmetry, the 
various sets of weights should have the same distribution, we arrive at the following saddle point equation: 

${h) = -2 ln(2.g) - In — ^ , for ft S [fti, ftz] . (22) 

n — hi 

This equation constitutes a standard Hilbert problem 21j whose solution is given by 

, r dz 2 ln(2g) + In 

H{h) = ^{h-h,){h-h,) j ^- , (23) 

2Tn [z - h)^/{z - hi){z - fta) 

with the contour 7 encircling counterclockwise the cut [fti,ft2]. The integral can be evaluated by inflating 
the contour and catching the pole at z = h and the logarithmic cut on [0, fti]. Taking the signature of the 
square root to be + to the right of the cut, and — to the left (and ±z respectively above and below the cut), 
and adding back the logarithmic part of the resolvent we find 



H{h) = In 



fti 



ft(fti + ft2) - 2ftift2 + 2\/ftIft^7(ft-fti)(ft-ft2) 



(24) 



On this expression we have to impose the condition -ff (ft) ^ ^ for large ft, which is a consequence of the 

normalization of the density function p(ft), and this fixes the parameters fti and ft2 to be hi — \^'^g2 and 

^2 — ]1ltg'^ ■ As expected we find a critical point at gc = 1/2, which is exactly the critical point of CDT [6i|. 
At the critical point fti = and ft2 = +00, i.e. the distribution has no saturated part and is excited all the 
way to ft = +00. 

To turn the above derivation into a proper solution it would of course be fundamental to prove the 
conjecture behind formula (|17p . but it is encouraging to see that this gives the right result for the critical 
point. 
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4 Conclusion 



Matrix model techniques proved to be a powerful tool in the study of 2D gravity before the introduction of 
CDTs. By combining ideas from multi-matrix models and dually weighted matrix models, we have been able 
to extend the technique to CDTs in a direct way. This has been achieved using the idea of "rigidity" : that 
constraints on local structures in a fatgraph can have global consequences. As we have seen, a restriction on 
colourings of vertices and faces serves to impose the CDT condition with all its consequences (for example, 
the restriction on global topology). It has been observed that the CDT model is intermediate between the 
unrestricted dynamical spacetime of the DT model, and the totally rigid regular lattice. It is interesting to 
review this observation in the light of the above results. 

The standard CDT model of the type discussed above rules out spatial topology change, and it is inter- 
esting to consider re-admitting topology change in a controlled way in the above model, as has been done in 
[12]. As mentioned in section [2] this can be done by allowing vertices connecting to four spacelike edges {e.g. 
by allowing vertices), or faces with four timelike edges (the dual of this is shown in figure [T]). The latter 
case would mean altering ([3]), setting Tr (C™) = NS2,m + NcSi^m for m> I, where c is a new "dual coupling" 
constant. Of course adding such extensions without control would bring us back to the universality class of 
the Euclidean DT, so we would instead have to scale correctly the new coupling in the continuum limit to 
control the topology change, in a similar fashion to what done in 12]. It would be interesting to compare 
the results of this approach with those of [TTl [T^l [13] and maybe establish a direct link. 

The matrix model above will be studied further in a forthcoming paper, using generalisations of the 
character expansion methods of |141 1151 116] . Hopefully, the power of matrix model techniques can then be 
applied to unsolved problems in 2D CDTs. Outstanding among these is the addition of matter, and the 
derivation of interesting results like those found for the "pure gravity" case. By comparison to previous work 
on DT matrix models, one could write down models for these cases without great difficulty. For example, 
one way of putting the Ising model on the CDT matrix model would be to double the number of matrices, 
replacing A with and A_ (and perhaps also B with and _B_) in an analogous way to the case of 
DTs. Application to the Potts model and others would be similar. 

The new perspective given here on CDTs is neither confined to 2D nor to CDT. In three and four 
dimensions also, analytic techniques that were of use in DTs [22] have previously been inapplicable to CDTs, 
due to the form of the causal condition. However, there are extensions of the claims of section [2] that can be 
conjectured to hold in higher dimensions. This leads to a reformulation of the causal condition. It would be 
most interesting to determine if this reformulation can overcome the barriers to analytic progress in CDTs in 
three and four dimensions. Investigations of this question by Daniele Oriti, Pedro Machado and the authors 
are currently underway. Furthermore these new way of imposing causality on the graphs might be of interest 
for other approaches to quantum gravity such as Group Field Theory |23j . 
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